First observations of a dome-shaped large-scale coronal EUV wave 
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ABSTRACT 



We present first observations of a dome-sliaped large-scale EUV coronal wave, 
recorded by the EUVI instrument onboard STEREO-B on January 17, 2010. 
The main arguments that the observed structure is the wave dome (and not 
the CME) are: a) the spherical form and sharpness of the dome's outer edge 
and the erupting CME loops observed inside the dome; b) the low-coronal wave 
signatures above the limb perfectly connecting to the on-disk signatures of the 
wave; c) the lateral extent of the expanding dome which is much larger than that 
of the coronal dimming; d) the associated high-frequency type II burst indicating 
shock formation low in the corona. The velocity of the upward expansion of the 
wave dome {v ~ 650 km s"^) is larger than that of the lateral expansion of the 
wave {v ~ 280 km s^^), indicating that the upward dome expansion is driven all 
the time, and thus depends on the CME speed, whereas in the lateral direction 
it is freely propagating after the CME lateral expansion stops. We also examine 
the evolution of the perturbation characteristics: First the perturbation profile 
steepens and the amplitude increases. Thereafter, the amplitude decreases with 
j,-2.5±o.3^ the width broadens, and the integral below the perturbation remains 
constant. Our findings are consistent with the spherical expansion and decay of 
a weakly shocked fast-mode MHD wave. 

Subject headings: Sun: corona — Sun: coronal mass ejections (CMEs) — Sun: flares 
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Introduction 



Large-scale propagating disturbances in the solar atmosphere occuring in association 



with flares and coronal mass ejections (CMEs 



filtergrams (JMoreton 



1960 



Athay &: Moreton 



hav e been first observed in chromospheric 



196ll ). These "Moreton waves" propagate 



with typical velocities in the order of 1000 km s~^, and have been interpreted as the 
ground-track of a dome-shaped MHD wave front propagating through the solar corona, 
wh ich conipresse s and pushes the chromospheric plasma downward when sweeping over 



it ([Uchida 



19681 ) ■ The Extreme-ultraviolet Imaging Telescope (EIT; 



Delaboudiniere et al 



19951 ) onboard S OHO for the fir s t time imaged such w ave-like disturbances in coronal 



emission lines ( jMoses et al. 



1997 



Thompson et al 



19981) 



EIT waves have been initially interpreted as the coronal counterparts of the 
chromo spheric Moreton waves as predicted in Uchida's fast-mode coronal MHD wave 



model (JThompson et al. 



19991 ). However, differences in morphology and propagation 



yelocities of Moreton waves and EIT wayes, w 



( iKlassen et al. 



2000 



Thompson fc Myers 



lich lie mostly in the range 200-400 km s ^ 



20091) led to severe doubts of this interpretation 



and alternative models were put forward. Some of them question if the phenomenon is a 
wave at all, and instead suggest that it is a signature of the large-scale coronal restructuring 
due to the erupting CME causing plasma compression or loca lized energy release (e.g. 



Delannee &: Aulanier 



1999 : 



Chen et al. 



2002; 



Attrill et al. 



the di fferen t models, we refer to t h e rece nt reviews by 



20071). For de t ailed discussions o f 



(120081), and 



Wills-Davev fc Attrilll f l2010h 



WarmuthI (120071 ) 



Vrsnak fc Oliver 



There seems to be some consensus that Moreton waves are indeed shock waves, as 
is suggested by their high propagation velocities and perturbation characteristics. It 
was shown that the amplitude and velocity of Moreton waves decrease, and the width 
of the wave pulse broadens as it propagates, consistent with a large-amplitude wave 
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or freely propagating shock wave that formed by steepening of a simple wave (e.g. 



Warmuth et al 



20041 ). However, the situation is unclear for large-scale waves observed 



in the corona. There is rather qualitative insight that the wave fronts not only become 



more diffuse but also broaden as they p ropagate (JThompson et al. 



2000 



Podladchikova fc Berghmans 



or dispersion of the wave (e.g. 



1999 



Klassen et al. 



20051). which may be du e to energy flux conservation 



different EIT waves, 



Wills-Davev fc Attrill 



20101 ). On a statistical basis of 



WarmuthI (l2010l ) concludes that for larger distances, the perturbation 



amplitudes tend to become smaller and the w idth larger 



evolution in high-cadence TRACE images by 



Wills-Davev 



oweyer, a case study of the wave 



20061) revealed that the width of 



the wave pulse remained constant during the propagation. 

Due to the low EIT cadence of ~12-15 min it was not possible to detect fast EIT 
waves and to study the evolution of the wave pulse characteristics, which is important in 



const raining the physical processes. The Extreme Ultraviolet Imager (EUVI; 



Howard et al 



20081 ) instr uments onboard the twin spacecraft of the Solar- Terrestrial Relations Observatory 



(STEREO; iKaiser et al.l 120081 ) overcome these limitations and offer several advantages 
for the study of large-scale coronal waves, in particular due to their high cadence, large 
field-of-view (FoV), high sensitivity, and the simultaneous observations from two vantage 
points. 

Since the launch of the STEREO satellites in October 2006, a variety of case 



studi e s of large-scale coronal waves have been carried out using EUVI data ( 



2008 



Veronig et a 



2008 



Gopalswamv et al 



2009 



Patsourakos &: Vourlidas 



ong et al. 



2009 



Patsoura 



cos et al. 



Ma et al. 



2009 



20091 



Attrill et a 



Zhukov et al. 



2009 



2009 



Cohen et al. 



Dai et al. 



20091 



Kienreich et al. 



2009; 



20101). The typical propagation velocities 



that were derived from EUVI waves (that all occurred during solar minimum conditions 
and were not accompanied by Moreton waves) are in the range 200-350 km s~^. For the 
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kinematical evolution, the results so far are rather inconclusive: some of the waves studied 
with EUVI showed evidence for deceleration during their propagation, consistent with the 



decay and deceleration of a 
coro na (e.g. 



[e.g. 



Long et al. 



Ma et al. 



2009 



2008 



arge amplitude wave to the fast-mode speed of the ambient 



Veronig et al. 



Kienreich et al 



20081 ). whereas others showed constant velocity 



20091 ). It is worth noting that many of the EUVI wave 



studies rev ealed a close a.ssociation of the wave and the erupting CME and its expanding 



flank s (e.g. 



Veronig et al. 



2008; 



Attrill et al. 



2009; 



Kienreich et al. 



2009 



Patsourakos et al 



20091 ). whereas the associated flares were all very weak. Stereoscopic EUVI studies revealed 
that the wave signal is typically confined to about 1-2 coronal scale heights above the 



chromosphere (jPatsourakos fc Vourlidas 
2OO9I). 



2009 



Patsourakos et al. 



2009 



Kienreich et al. 



In this letter, we present the first observations of the full dome of the wave observed 
in EUV, which is consistent with the three-dimensional expansion of a coronal shock front. 



We note that part of a wave dome had been observed in soft X-rays by Yohkoh/SX' 



fast coronal wave that occurred in association with a Moreton wave (iNarukage et al. 



for a 



2004h . 



We also study the evolution of the wave pulse characteristics (velocity, amplitude, width) 
in high-cadence EUV imaging over a propagation distance of more than IRq, and discuss 
the implications for the physical processes involved. 

2. Data 

The EUVI instrume nt is part of the Su n Earth Connection Coronal and Heliospheric 



Investigation (SECCHI; 



Howard et al 



20081 ) instrument suite onboard the STEREO- 
A(head) and STEREO-B(ehind) spacecraft. EUVI observes the chromosphere and low 
corona in four spectral channels (He ll 304 A: T ~ 0.07 MK, Fe ix 171 A: T ~ 1 MK, 
Fe XII 195 A: T ~ 1.5 MK, Fe XV 284 A: T ~ 2.25 MK) out to 1.7 Rq with a pixel limited 
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spatial resolution of 1.6"/pixel ( IWuelser et al 



20041 ) ■ On 2010 January 17, STEREO-B was 



69.2° behind Earth on its orbit around the Sun, observing a large-scale coronal wave in 
its Eastern hemisphere. The EUVI-B imaging cadence was 2.5 min in the 171 A, 5 min 
in the 195 A, 2.5-5 min in the 284 A, and 5 min in the 304 A passband. The EUVI data 
were reduced using the secchi_prep routine available within Solarsoft, and corrected for 
solar differential rotation before we derived base and running difference and ratio images, 
respectively. We also show white light observations from the CORl-B coronagraph, which 
has a FoV of 1.5-4 R©, and observed with a cadence of 5 min during the event under study. 

3. Results 

Figure [1] shows the evolution of the wave in EUVI 195 A direct and running difference 
images, where we subtracted from each frame the frame taken 5 min before (see also the 
accompanying movie no. 1). The wave is best observed in the EUVI 195 A filter, which 
has a broad temperature response peaking at ~1.5 MK, but can be observed in all four 
wavelengths. In the 195 A images at 03:56 and 04:01 UT the full dome of the wave is clearly 
observed, even in the direct images. The image sequence also reveals that the on-disk 
signatures of the wave perfectly connect to the wave dome observed above the limb. The 
sharp and very regular edges of the dome further suggest that we really observe the shock 
front of the wave. In EUVI 171 A images, erupting CME loops are observed behind/inside 
the dome (cf. Fig. [2]). 

In Fig. [2], we plot difference images taken almost simultaneously (between 03:56 
and 03:57 UT) in all four EUVI channels. The wave dome can be identified in all four 
wavelengths, which indicates that structures at different temperatures are disturbed by the 
wave, covering at least the temperature range 1.0-2.3 MK. The fact that we can observe 
the wave dome also in the EUVI 304 A filter actually indicates that Si xi is significantly 
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cont ributing to 



;he 304 A emission in addition to the 10^ K emission due to He ll lines (see 



also 



Long et al 



2008 



Patsourakos fc Vourlidas 



20091) 



We also note that HIRAS (Hiraiso Radio Spectrograph) reported an associated 
high-frequency type II burst drifting from ~310 MHz to ~80 MHz during ~03:51-03:58 UT. 
The wave center was derived at a meridional distance of 57° for STEREO-B (cf. Fig. ^, 
which implies that it was located 36° behind the Eastern solar limb for an Earth-based 
vantage point, corresponding to an occultation height of about 0.23 i?© ~ 160 Mm. Since 
the radio source was behind the solar limb when looking from the Earth, the observed 
emission has to be at the harmonic of the plasma frequency. Applying two- to ten-fold Saito 
density models this corresponds to the height range 0.11-0.35 R©, suggesting that the shock 
occurred relatively low in the solar corona. These shock formation heights are consistent 
with the heights of the wave dome observed in EUV (cf. Fig. H]). 

In the sixth panel of Fig. [H the outer contours of the coronal dimming region as 
identified in 195 A base ratio images (05:01 UT/03:36UT) are overlaid. At the time, at 
which we extract the contours, the dimming was maximally developed and darkest; the 
contour lines plotted in Fig. [1] are at 80% of the pre-event intensity. It is evident that the 
North-South extent of the dimming region, which outlines the lateral extent of the CME 
structure low in the corona (CME flanks), is significantly smaller than that of the expanding 
wave dome. If the dome corresponded to the CME body, its extent should not exceed that 
of the coronal dimming. In Fig. |3] (panels b and c) we show composits of EUVI 195 A 
and CORl images. The wave dome observed in EUVI images smoothly extends to the 
white- light structure observed in CORl. This suggests that the outer edge observed in the 
white-light coronagraphic images correspond rather to the coronal shock ahead of the CME 



than to the CME leading edge itself (see also 



Vourlidas et al 



20031) 



Figure |3^ shows an EUVI-B 195 A difference image together with all wavefronts 



determined in the 195 A passband and t he center obtained f rom circular fits to the earhest 
wavefronts in the 3D spherical plane (cf. 



Veronig et al. 



20061 ) ■ For each wavefront visually 



identified in the 195, 171, 284 and 304 A difference images, we determined the mean 
distance from the derived center along the spherical solar surface. The top panel in Fig. H] 
shows the resulting wave kinematics. The velocity of the wave obtained from the linear fit 
to the kinematical curve is f ~ 283 ± 27 km s^^, and remains constant over the propagation 
distance up to 950 Mm. In the same panel we also plot the kinematics of the wave dome 
followed along its main propagation direction as observed in EUVI and CORl. The distance 
of the wave dome is measured against the plane-of-sky, and for the starting point we use 
the center derived for the on-disk EUVI wave measurements. The plot shows that at the 
beginning the distance of the on-disk wave and the height of the dome are roughly in 
agreement, but the upward movement of the wave dome is much faster {v ~ 650 km s~^) 
than the lateral expansion of the wave observed on the solar disk {v ~ 280 km s^^). 

In Fig. [5] we plot the evolution of the intensity amplitudes of the wave, so-called 
"perturbation profiles" , determined from 195 A ratio images, where we divided each frame 
by the frame recorded 10 min before. We calculated the intensity profiles over a 60° sector 
on the solar sphere (indicated in Fig. [3^), where the signal of the wave is strongest, by 
averaging the inten sities of all pixels in "rings" of increasing radius arou nd the wave origin 



shown in Fig. |3] (cf. 



Podladchikova fc Berghmans 



2005 



Muhr et al 



20101 ). Base ratio images 



would in principle be better suited to study the perturbation profiles. However, these are 
affected by changes in the quiet Sun as well as by brightenings induced by the wa ve front 



passage that are only slowly fading ( "stationary brightenings" ; lAttrill et al.l 



20071 ). As a 



result of several tests, we use 10-min running ratio images since they well represent the 
propagating wave profile and ensure that the peak of the wave amplitude is not cut. 



The propagation of the wave is well refiected in the perturbation profile evolution 
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plotted in Fig. [5] (see also the accompanying movie no. 2). We also observe steepening and 
amplitude increase in the profiles 03:56 to 04:01 UT, where the highest amplitude of 1.45 
is reached, corresponding to an enhancement of 45% above the pre-event level. Assuming 
that the intensity enhancement is primarily due to plasma compression rather than due to 
temperature changes, which is somewhat justified by the observations of the wave over a 
broad temperature range in the four EUVI channels, this intensity amplitude of ///q = 1.45 
in an optically thin emission line corresponds to a density ratio n/riQ oc {I/IoY^"^ ~ 1.2. 
The perturbation amplitudes in the other EUVI channels are smaller than in 195 A but 
clearly recognized; see the peak amplitudes in 171 A and 284 A plotted in the second panel 
of Fig. [5l We also note that at the time of maximum amplitude, the outer edge of the wave 
front is steepest. In the subsequent evolution, the amplitude and the steepness of the wave 
front decrease until it can no longer be followed in the profiles at about 04:36 UT. 

In the middle panel of Fig. H] we plot the evolution of the intensity amplitude, 
which decreases continuously after the peak at 04:01 UT. We fitted the profile 
amplitude A = {I — /o)//o as a function of the propagation distance d with a power-law 
A = a ■ d'^, giving b = —2.5 ± 0.3. We also measured the width of the wave and followed its 
evolution. Due to the problematics of the stationary brightenings behind the wave front 
(which masks the trailing part), we extract only the width of the frontal part, i.e. from 
the profile peak to the outer front. We derived the full width of the frontal part (defined 
from the peak down to 2% enhancement above the pre-event level), as well as the frontal 
width at half maximum. These are plotted in the bottom panel of Fig. H] together with the 
integral of the perturbation profile (as shown in Fig. E]) over the frontal part of the wave. 
The width of the wave pulse increases during its evolution by a factor of 3-4, whereas the 
integral remains basically constant (changes are <25%). 
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4. Discussion and Conclusions 

There are four main arguments supporting tlie conclusion tliat tlie observed EUV 
structure is tlie wave dome and not the CME body: a) The dome appears spherically 
three-dimensional with a sharp outer edge; inside the dome erupting CME loops are 
observed, b) The wave signatures observed in the low corona above the limb perfectly 
connect to the wave signatures observed against the solar disk, c) The lateral extent of 
the expanding dome is much larger than that of the coronal dimming, d) The event was 
associated with a high-frequency type II burst indicating shock formation low in the corona, 
which is consistent with the observed EUV dome height. 

The velocity derived for the upward expansion of the wave dome {v ~ 650 km s~^) is 
larger than that of the lateral expansion of the wave observed on-disk {v ~ 280 km s^^). 
There are at least two alternative explanations for this. In the freely propagating phase, 
the upward-lateral velocity difference can be due to differences of the fast magnetosonic 
velocities of the active region (AR) and the ambient quiet corona. The tip of the wave 
dome lies above the AR, and the m agnetosonic speed is muc h higher above ARs than at 



low heights in its surroundings (e.g. 



Warmuth fc Mann 



20051 ). An alternative explanation 



is that the upward dome expansion is driven all the time, and thus depends on the 
CME speed, whereas in the lateral direction the wave is freely propagating as soon as 
the lateral expansion of the CME flanks has stopped, and its velocity is determined by 
the characteristic speed of the medium. The second interpretation is supported by the 
evolution of the dimming region in the low corona, which is fast expanding up to about 
04:01-04:06 UT. Thereafter, the coronal dimming still gets darker until about 05:00 UT, 
indicating ongoing mass depletion, but with only little (or no) further expansion after 
04:06 UT. 

The high-cadence EUV observations together with the distinct wave signal observed 
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in the low corona against the solar disk allowed us to study in detail the evolution of 
the perturbation characteristics during the coronal wave propagation. We find that the 
amplitude of the perturbation first increases and the perturbation profile steepens within 
the first 5 niin of the event (i.e. until 04:01 UT), thereafter the amplitude A decreases, 
following a power-law of the form A oc r~^'^''=°'^. The width of the perturbation profile 
broadens during its evolution by a factor of 3-4. This broadening is observed for both the 
full width as well as for the width at half maximum. We stress that we only measured the 
frontal width of the wave pulse due to the uncertainties in the trailing part (mostly due to 
stationary brightenings induced by the wave passage). Since the perturbation profile is not 
necessarily symmetric, the evolution of the total width of the wave pulse may be different, 
but we do not expect that it would alter the main outcome of broadening, which is a quite 
pronounced effect. The integral below the frontal part of the perturbation profile basically 
remains constant. These characteristics of the wave pulse evolution are all consistent with 
the spherical expansion and decay of a large-amplitude fast-mode MHD wave. 

The velocity of the lateral wave expansion {v ~ 280 km/s) lies in the range of the 



fast magneto sonic speed in t 



discussion in 



le qu iet solar corona during solar minimum condition (see 



Kienreich et al.l 



20091 ). From the observed EUV wave intensities at the highest 



amplitude reached at ~04:01 UT, we obtained a rough estimate of the peak density jump 



at the wave front, n/nQ ~ 1.2. 



number of Mfms ~ 1.15 ((Priest 



lis co rresponds to a perpendicular fast magnetosonic wave 



19821 ). Thereafter, as the wave decays to n/uQ = 1.04 
(Fig. [5l last profile at 04:36 UT), it approaches a hnear regime since Mfms ~ 1-03. These 
numbers indicate the evolution of a weak shock, and may explain why we do not observe a 
significant deceleration of the free-propagation of the lateral wave expansion in association 
with the amplitude decay. 

AMV, NM, and IWK acknowledge the Austrian Science Fund (FWF): P20867-N16. 
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Fig. 1. — Sequence of STEREO-B EUVI 195 A images showing the early wave evolution on 
the disk and above the limb (top panels: direct images, middle and bottom images: 5-min 
running difference images). [See also accompanying movie no. 1.] Note the dome shape in 
the images at 03:56 and 04:01 UT. In the sixth panel, we overlay the outer contours of the 
coronal dimming. 
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Fig. 2. — 5-min running difference images of tlie wave dome in all four EUVI-B spectral 
channels. Arrows outline the wave dome; crosses indicate the erupting CME loops inside the 
dome. 
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Fig. 3. — a) Median-filtered EUVI-B 195 A 5-min running difference image at 04:06 UT. The 
yellow curves indicate the outer edges of the wavefronts visually identified in the EUVI-B 
195 A filtergrams. The blue curves outline the wave propagation sector in which the intensity 
profiles plotted in Fig.Elwere calculated, starting from the determined wave center (indicated 
by a yellow star), b), c) Composits of EUVI-B and CORl-B running difference images (note 
the larger FoV compared to panel a). The black circle indicates the inner border of the 
CORl FoV at 1.5 Rq. 
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Fig. 4. — Top: Kinematics of the wavefronts observed on the solar disk in all four EUVI- 
B channels together with the upward motion of the wave dome measured in EUVI-B and 
CORl-B. Middle: Evolution of the perturbation amplitude determined from the EUVI-B 
195 A intensity profiles shown in Fig. O Bottom: Evolution of the width (at half maximum: 
diamonds, full: triangles) and the integral (stars) of the frontal part of the perturbation 
profiles. 
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Fig. 5. — 195 A intensity enhancements ("perturbation profiles") of tlie propagating wave 
derived over a 60° sector wliere tlie wave is strongest (cf. Fig. [3^). Note tliat all values 
behind the propagating wave front which are smaller than one (as a result of the running 
ratio procedure) are set to one. In the second panel, we plot also the peak intensity profiles 
derived from the EUVl-B 171 and 284 A images. [In the accompanying movie no. 2, we 
show the wave evolution in the EUVI-B 195 A 10-min running ratio images together with 
the corresponding perturbation profiles.] 



